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Introduction: Recent studies have shown that 

Martian magmas had wide range of oxygen fugacities 
(f0 2 ) and that this variation is correlated with the 
variation of La/Yb ratio and isotopic characteristics of 
the Martian basalts, shergottite meteorites [1,2]. The 
origin of this correlation must have important 
information about mantle sources and Martian 
evolution. In order to understand this correlation, it is 
necessary to know accurate value of oxidation state of 
other Martian meteorite groups. Nakhlites, cumulate 
clinopyroxenites, are another major group of Martian 
meteorites and have distinctly different trace element 
and isotopic characteristics from shergottites [3]. Thus, 
estimates of oxidation state of nakhlites will give us 
important insight into the mantle source in general. 

Several workers have estimated oxidation state of 
nakhlites by using Fe-Ti oxide equilibrium [e.g. 4,5]. 
However, Fe-Ti oxides may not record the oxidation 
state of the parent melt of the nakhlite because it is a 
late-stage mineral. Furthermore, there is no 
comprehensive study which analyzed all nakhlite 
samples at the same time. Therefore, in this study (1) 
we reduced the uncertainty of the estimate using the 
same electron microprobe and the same standards 
under the same condition for Fe-Ti oxide in 6 nakhlites 
and (2) we also performed crystallization experiments 
to measure partition coefficients of Eu into pyroxene 
in the nakhlite system in order to estimate f0 2 when 
the pyroxene core formed (i.e. Eu oxybarometer [e.g. 
2 , 6 ]). 



Fig. 1 BSE images of Fe-Ti oxides in Nakhla (left) and 
in MIL03346 (right). The bright areas are magnetite and 
light gray areas are ilmenite, which are surrounded by augite 
or mesostasis. We can clearly see wider ilmenite lamellae 
(~20 pm across) in Nakhla compared to much thinner 
lamellae (< 1 jam across) in MIL03346. 


Method: Fe-Ti oxides (magnetite-ilmenite) in 6 

nakhlites (MIL03346, NWA817, Y000593, Nakhla, 
Govemador Valadares, NWA998) were measured. We 
calculated their f0 2 using two different 
geothermometer/oxygen models, Ghiorso and Sack [7] 
and the Ca-QUILF [8] models that differ in elements 
accounted for as well as in the solution models used. 

We calibrated the Eu oxybarometer for the nakhlite 
using the synthetic NJ4 composition, which we believe 
is the closest to the Nakhla parent melt [9]. We doped 
this NJ4 composition with percent level REEs... Then, 
we put this starting material in the furnace at 1300 °C 
for 48 hours to homogenize it, quenched and put back 
below the liquidus (1150 ~ 1160°C) and grew 
pyroxenes isothermally or during slow cooling. Runs 
were performed under the three oxygen fugacities 
between the Iron-Wustite (IW) buffer and the Quartz- 
Fayalite-Magnetite (QFM) buffer (~ 3.5 log units 
above IW). We analyzed the resulting samples with 
the JEOL JXA8900L electron microprobe at 
University of Tokyo by using well-characterized REE 
standards and calculated the D(Eu/Gd). 


Table 1 . Result of f02 estimates for the nakhlite by the Fe-Ti method 



TOC) 

logfOZ (atm) 

log fOj 
(A QFM) 

Gliiorso & Sack [7] 

MIL-0334 6 

n/a 

n/a 

n/a 

N WAS 17 

618-719 

-16.92 to -20.61 

033 to -0.72 

Y000593 

610-658 j 

1-19.69 to -21 34 

-133 to -1.19 

G.V. 

655-692 

-18.50 to -20.04 

-034 to -129 

NnkliJh 

648-720 j 

| -1730 to -19.93 

-039 to -139 

NWA99S 

674-680 

-1831 to -18.71 

-026 to -032 

Ca-QUILF [SI 

MIU33346 

n/a 

n/a 

n/a 


667-724 j 

! -16.46 to -1834 

-037 to -038 


650-680 

-1838 to -1937 

! -130 to -1.14 

G.V. 

649-688 

j -1820 to -19.85 j 

j -031 to -1.40 

Nakhla 

671-712 

-17.70 to -1831 

-038 to -1.10 

NWAPP8 

682 - 700 

i-1730 to -1839 | 

| -027 to -038 


Results and Discussion: Ilmenite lamellae in 

MIL03346 were too thin (<1 jam) for analysis by 
microprobe (Fig. 1). The width of ilmenite lamellae in 
nakhlites studied is consistent with the cooling model 
proposed for the nakhlite igneous body [10]. We 
checked whether magnetite and ilmenite are in 


equilibrium by using the Mg/Mn partitioning test from 
[11]. All Fe-Ti oxides we chose have cation sums 
within error of 3.00 and are therefore in stoichiometry. 
The f0 2 results by the Fe-Ti oxide method are 
summarized in Table 1. The absolute error on the f0 2 
estimates is ±0.5 log unit at most as reported in [12]. 
The estimated values are quite similar in both models. 
Although the f0 2 estimates for NWA817 and 
NWA998 are slightly higher than those for the other 
nakhlites, there is generally good agreement among the 
f0 2 estimates and equilibrated temperatures of the 5 
nakhlites within error. This result suggests that Fe-Ti 
oxides in all 5 nakhlites were equilibrated at 650 ~ 700 
°C under QFM-0.5 ~ QFM-1. This oxidation state is 
consistent with the previous estimates using Fe-Ti 
oxide equilibrium by several authors [e.g. 4,5]. 



log f02 QFM) 


Fig. 2 Calibration of D(Eu)/D(Gd) in synthetic pyroxene 
versus oxygen fugacity (in log unit relative to QFM buffer). 


lower by 1-2 log unit than those estimated from the 
Fe-Ti oxides [2,13]. The correlation between the 
whole rock La/Yb ratio and f0 2 seen in shergottites 
seems not to be applicable to the nakhlite group (Fig. 
3), indicating that the nakhlite source had a very 
different evolutional history than the mantle from 
which the shergottites were derived. 


Table 2. Pyroxene D(Eu)/D(Gd) 

Expt 

f0 2 | 

(A QFM) j 

D(Px) 

Nakltla melt 

817 

-33 ! 

037 ±0D1 


818 

-2D ! 

0j66 ±0D1 

0.70 ±0D4 

819 

0 ! 

0.74 ±0D1 




Fig. 2 shows our calibration of D(Eu/Gd) in 
synthetic pyroxene vs. oxygen fugacity for our 
nakhlite experiments. This calibration curve has a 
similar trend to that for the shergottite calibration [13] 
and that predicted by theory [6]. Unlike the shergottite, 
the bulk composition of the nakhlite dose not represent 
its parent melt composition. Thus, in order to calculate 
D(Eu)/D(Gd) of nakhlite melt, we used the partition 
coefficients of REE in the synthetic Nakhla augite/melt 
from the experimental study [14] as well as the 
published REE data of Nakhla augite core [15] and 
Nakhla whole rock [16]. Table 2 shows this calculated 
D(Eu)/D(Gd) value of Nakhla melt and our 
experimental D value under 3 oxygen fugacities. The 
f0 2 for Nakhla was estimated by comparing these two 
D values. Our study suggests that Nakhla augite 
crystallized at 1150-1160 °C under QFM-1 (±1.0) 
buffer (Table 2 and Fig.2). 

These two independent results show us that the 
magmatic oxidation state for the nakhlite did not 
change throughout their crystallization. This presents a 
contrast to the case of the shergottite, in which the f0 2 
values using the Eu oxybarometer are systematically 


Fig. 3 Whole rock La/Yb ratio versus oxygen fugacity for 
the martian meteorites. The f0 2 data are from [1,17,18] and 
La/Yb data are from [3,17,19,20]. 
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